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The elliptic flow {V 2 ) at -^snn = 11.5, 39, and 200 GeV and triangular flow (113) at .^snn = 200 
GeV of identified particles , Kg,p,p,(j>, A and A) from 0%-80% central Au+Au collisions 

are analyzed using a multiphase transport (AMPT) model. It is shown that the experimental results 
from the rj-suh event plane method can be reproduced with a parton scattering cross section between 
1.5 and 3 mb. We also studied the differential and integrated V 2 /V 3 ratios and conclude that they 
are anti-correlated with the parton scattering cross section. 

PACS numbers: 25.75.-q, 25.75.Ld 


I. INTRODUCTION 

One of the main goals of heavy-ion collision experi¬ 
ments at the Relativistic Heavy Ion Collider (RHIC) is 
to study the properties and formation of the Quark Gluon 
Plasma (QGP). The study of the azimuthal anisotropy 
of emitted particles, based on Fourier decomposition, is 
considered to be one of the most important tools to in¬ 
vestigate the hot and dense matter created in heavy-ion 
collisions [T]. Several interesting observations of the sec¬ 
ond harmonic Fourier coefficient V 2 , also called elliptic 
flow, have been reported during the past decade by us¬ 
ing the data from the top RHIC heavy-ion energy of 
Y^snn = 200 GeV [5H1]. At low transverse momenta 
(pt < 2.0 GeV/c), a mass ordering of V 2 {pt) was ob¬ 
served mm, which can be understood within a hydro- 
dynamic framework [7]. At intermediate transverse mo¬ 
menta (2.0 < Pt < 6.0 GeV/c), a Number-of-Gonstituent 
Quark (NGQ) scaling [8] was observed [IHIl]. The NGQ 
scaling was used to conclude that the relevant degrees of 
freedom in the created matter are quarks and gluons. In 
addition to measurements at top RHIG energies, a Beam 
Energy Scan (BES) program has been carried out in the 
years 2010, 2011 and 2014 at to study the structure of the 
QCD phase diagram. Au-|-Au collisions were recorded at 
yisWN = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV. 
The mass ordering in the low px region and NGQ scaling 
at intermediate transverse momenta were also observed 
BES energies [T?lfli] . 

In the recent years, the third harmonic Fourier coeffi¬ 
cient U3, called triangular flow, has attracted more atten¬ 
tion and was studied intensively nsjus]. It is thought 
to be related to the near-side ridge structure observed 
in two particle correlation measurements |15j . The tri¬ 
angular flow seems to be more sensitive to the viscous 
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effects than V 2 m, but less sensitive to the collision cen¬ 
trality m- Thus, U3 is a natural choice to study initial 
fluctuations and viscosity effects. Similar flow patterns as 
for V 2 were recently observed in measurements HHHO]. 

Viscous relativistic hydrodynamic calculations unm] 
suggest a combined analysis of elliptic and triangular 
flow data to put tight constraints on the initial states 
of nucleus-nucleus collisions. Other hydrodynamic model 
calculations [23] predict a stronger damping of the higher 
harmonic flow coefficients (n > 2) relative to V 2 at high 
transverse momenta. A constant ratio of U2 (pt)/i' 3 (pt) 
at high Pt is predicted, where high means px /tot > Vmax 
and Vmax is the highest flow velocity. Thus, by combin¬ 
ing elliptic and triangular flow, more information about 
the initial state and the influence of the viscous effect can 
be collected. 

In this paper, the differential and integrated ratio of 
V 2 /V 3 of identified particles {tt^, , Kg,p,p,(j), A and 

A) are studied with the AMPT model. This paper is 
organized as follows. Section [IT] gives a brief introduction 
to the 77-sub event plane method and the AMPT model. 
In section HI A the collective flow calculations {v 2 and 
V 3 ) from the AMPT model are presented. The differential 
and integrated V 2 /V 3 ratios are discussed in Sec. HIB and 
sec. inrcj A summary is given in Sec. |TV1 


II. 77-SUB EVENT PLANE METHOD AND THE 
AMPT MODEL 

A. 77-sub Event Method 

The event plane method is one of the most widely 
used methods to analyze anisotropic flow in heavy-ion 
collisions [23]. In non-central Au-|-Au collision, the over¬ 
lap region has an almond shape in the coordinate space. 
As the system evolves, the pressure gradient pushes the 
anisotropy from coordinate space to momentum space. 
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therefore the produced particles have an anisotropic dis¬ 
tribution in momentum space. The azimuthal distribu¬ 
tion of the produced particles can be written as [IS] : 


E 


d^N 

dp^ 


1 d^N 

27r pT^dp^dy 


(1 + X! cos[n((?i-«'„)]), 

n—1 


( 1 ) 


where ij) is the azimuthal angle of a particle, is the 
observed n-th harmonic flow., and is the n-th har¬ 
monic event plane angle reconstructed by the produced 
particles, defined as: 




n 


1.1 Wj sm{n(j),) 
n Ei cos(u(/)i) 


( 2 ) 


The sum goes over all charged particles in the event, Wi is 
a weight applied to optimize the event plane resolution, 
see [23] for details. It is worth to note that the range of 
n-th harmonic event plane angle is 0 ^ 'I'n < 27r/n. 

In this method the n-th observed harmonic flow can 
be corrected for statistical effects by the n-th harmonic 
event plane resolution: 

(3) 


with 


Rn = (cos[n(5'„ - . (4) 

Here (...) denotes the average over all particles in all 
events and 'I'„u is the n-th real event plane angle pS] . 
The real event plane means the participant plane 
which can not be achieved due to finite multiplicity of an 
event. 

To reduce the short range ’’non-flow” effects such as 
HBT correlations, a Ty-sub event plane method was in¬ 
troduced [14]. The Ty-sub event plane method divides 
the event into two independent sub-events with different 
pseudo rapidity (ry) ranges. An additional gap of 0.1 in 
pseudo rapidity is added to further avoid short range cor¬ 
relations. The procedure results in a positive and neg¬ 
ative ry-sub event plane. The n-th observed harmonic 
flow can now be calculated in respect to the ?y-sub event 
planes. As a consequence, the ry-sub event plane resolu¬ 
tion [24] is used instead of full event plane resolution: 

= y/(cos[rr(4-“-vI/fc)]), (5) 

where the 4/“ and the are the n-th harmonic event 
plane angles for the positive and negative ry-sub ranges. 


AMPT (AMPT Def in following figures) and string melt¬ 
ing AMPT (AMPT Str in following figures). The initial 
conditions are generated by the HIJING (Heavy Ion Jet 
Interaction Generator) model p5II3Ij . HIJING includes 
only two body nucleon-nucleon interactions and gener¬ 
ates mini jets and excited strings through hard processes 
and soft processes separately. 

Excited strings are treated differently in the two 
AMPT versions. In default AMPT, excited strings are 
combined to hadrons according to the Lund string frag¬ 
mentation model they further go through a hadronic in¬ 
teraction stage [55]. In string melting AMPT, excited 
strings first convert to partons, i.e. melting, and then 
go through a partonic interaction stage with original 
soft partons. The partonic interactions for both default 
AMPT and string melting AMPT are described by ZPG 
(Zhang’s Parton Cascade) model [35]. In the final stage 
of the ZPG model, partons in default AMPT are recom¬ 
bined with parent strings and hadronize via the Lund 
string fragmentation model. However, in string melt¬ 
ing AMPT, the hadronization of partons is described by 
a coalescence model. After hadronization, the hadronic 
interactions are modelled by ART (A Relativistic Trans¬ 
port) model [531153] . 

In previous studies it was found that a large parton 
scattering cross section (approximately 6 mb or 10 mb) is 
needed to reproduce the observed elliptic flow of charged 
hadrons [351[3n]. In these studies, the elliptic flow was 
calculated with respect to the reaction plane (determined 
by the beam axis and impact parameter direction), which 
is 0 in the AMPT model, and this might underestimate 
the elliptic flow observed in the experiment data, which 
measured relative to the event plane but not the reaction 
plane [24] . As shown in Ref. m. both the charged parti¬ 
cle multiplicity and elliptic flow can be reproduced with a 
smaller parton scattering cross section (1.5 mb) by using 
the default values for the parameters in the Lund string 
fragmentation function. 

In this paper, approximately 20 million events for 0%- 
80% central Au-|-Au collisions at ^snn = 11-5, 39, and 
200 GeV with default AMPT (vl.25) and string melting 
AMPT (v2.25) were generated. Three different parton 
scattering cross sections (1.5 mb, 3 mb, and 6 mb) are 
used in the string melting version of AMPT. The param¬ 
eters in the Lund string fragmentation function are set 
to the default values for AMPT with 1.5 mb parton scat¬ 
tering cross section (the values can be found at Table I 
in Ref. |37j ). Furthermore the ry-sub event plane method 
is used to calculate elliptic and triangular flow. 


B. AMPT Model 

AMPT is a transport model which consists of four 
main components: the initial conditions, partonic in¬ 
teractions, conversion from partonic to hadronic mat¬ 
ter, and hadronic interactions |28j . It has two ver¬ 
sions to deal with different scenarios which are default 


III. RESULTS 

Identical analysis steps are used for simulation and 
data collected from the experiments [T4l[5D]. Finial state 
charged particles with 0.2 < pT < 2.0 GeV/c and |ry| < 
1.0, are used to reconstruct the event plane. The p range 
of the two ry—sub events are -1.0 < ry < -0.05 and 0.05 
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FIG. 1. (Color online) Second harmonic (top panel) and third harmonic (bottom panel) event plane resolution from the AMPT 
string melting model with 1.5 mb (AMPT Str (1.5 mb)), 3 mb (AMPT Str (3 mb)) and with 6 mb (AMPT Str (6 mb)) parton 
scattering cross section from 0%-80% central Au+Au collisions at .^snn = 11.5, 39 and 200 GeV. The analysis method is r/suh 
event plane method, and the event plane is reconstructed by charged particles, see text for details. 


< rj < 1.0, respectively. The elliptic and triangular flow 
are calculated within - 1.0 < 77 < 1 . 0 , and the auto¬ 
correlation is naturally avoided by the Ty-sub event plane 
method. For all the plots in this section, only statistical 
errors are shown. 

A. Elliptic Flow and Triangular Flow at RHIC 

Figure [^presents the second and third harmonic event 
plane resolutions from the AMPT string melting model 
for parton scattering cross sections of 1.5 mb, 3 mb and 
6 mb from 0%-80% central Au-I-Au collisions at ^snn = 
11.5, 39 and 200 GeV by using the ry-sub event plane 
method. The event plane resolution shows a clear parton 
scattering cross section and energy dependence: both sec¬ 
ond and third harmonic event plane resolution increases 
with increasing parton scattering cross section and en¬ 
ergy. 

The third harmonic event plane resolution is about two 
times smaller and peaks more central than the second 
harmonic event plane resolution, which is comparable to 
the one from experiment. For = 11.5 and 39 GeV, 

U 3 can only be calculated up to 0%-60% (0%-70%) for 
some of the parton scattering cross sections, since we only 
have event plane resolution up to these centrality. 

Figure E Fig- i and Fig. show the comparison be¬ 
tween data from experiment and from the AMPT model 
for V 2 {pt) of identified particles ,p,p and (j)) 


at ^snn = 11-5, 39 and 200 GeV from 0%-80% cen¬ 
tral Au-|-Au collisions. The data from experiment are 
taken from Ref. [H EH]- For all three energies it is 
observed that the elliptic flow increases with increasing 
parton scattering cross section. The elliptic flow from 
default AMPT is lower than the corresponding results 
from experiment, which indicates that the parton scat¬ 
tering process is important to produce a large elliptic 
flow as observed in experiment data. Previous studies 
show that a 6-10 mb parton scattering cross section is 
needed to reproduce the elliptic flow observed in experi¬ 
ment data [SS] [H] . In those studies the elliptic flow was 
calculated relative to the reaction plane. 

By using the ry-sub event plane method most of the 
experiment data points at low transverse momenta are 
between AMPT default and AMPT string melting with 
a 3 mb parton scattering cross section. At intermedi¬ 
ate pt some particle species, e.g. 7 r=*=, are systematically 
above AMPT string melting with a 3 mb parton scatter¬ 
ing cross section. That means previous studies have un¬ 
derestimated the elliptic flow by using the reaction plane 
instead of the event plane. Latter one fluctuates event- 
by-event around the reaction plane. The collective flow 
measured in experiment are not mean values, but closer 
to a root-mean-sqaure [53| , therefore, the event-by-event 
fluctuations did not cancel each other but give a posi¬ 
tive contribution to the measured collective flow, which 
makes the collective flow relative to the event plane is 
always larger compare to the collective flow relative to 
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FIG. 2. (Color online) V 2 {pt) results at ^snn = 11-5 
GeV from 0%-80% central Au+Au collisions for identi¬ 
fied particles (tt^ , A'* ,p,p and (j)) in comparison to STAR 
data. 



FIG. 3. (Color online) V 2 {pt) results at ^snn = 39 
GeV from 0%-80% central Au-fAu collisions for identi¬ 
fied particles (tt^ , , p, p and (p) in comparison to STAR 

data. 
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FIG. 4. (Color online) V 2 (pt) results at y^SNN = 200 FIG. 5. (Color online) V 3 {pt) results at y^SNN = 200 GeV 

GeV from 0%-80% central Au-|-Au collisions for identified from 0%-80% central Au-|-Au collisions for identified particles 
particles ,p,p and (p) in comparison to RHIC data. ,p,p and (p) in comparison to STAR data. 

The open symbols in the panels of ,p and p are from 

PHENIX |38| and the panel of cp are from STAR |10| and the 
lines are from AMPT calculation. 

the reaction plane. 

Figure shows the comparison between experiment 
data m and AMPT model calculations for V 3 {pt) at 
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y'SNN = 200 GeV from 0%-80% central Au+Au collisions 
of identified particles (tt^ , ,P,P and (j)). A similar par- 
ton scattering cross section dependence is observed as for 
the elliptic flow, but here the dependence is stronger: V 3 
is about 40% lower than V 2 with 1.5 mb parton scattering 
cross section but comparable with V 2 with a 6 mb parton 
scattering cross section. The driving force behind V 2 is 
the almond shape of initial nuclei overlap region, there¬ 
fore, V 2 has a relative large value compare to V 3 which is 
mainly generated by initial state fluctuations [H]. This 
effect can be observed by comparing the V 2 and V 3 val¬ 
ues from AMPT default model calculations in Fig. and 
Fig.i The higher parton scattering cross section makes 
the transportation from the initial coordinate space to 
the final momentum space more efficient, which increases 
the values for both V 2 and V 3 . The relative increase of 
V 3 is larger than for U 2 , since the effects of the parton 
scattering is more important for V 3 . The triangular and 
elliptic flow from experiment can be well reproduced by 
AMPT default and AMPT string melting with a 1.5 mb 
parton scattering cross section, see Fig. and Fig. 

It is worth to be noticed here that we did not find out 
a single parton scattering cross section to describe all 
particle species. This is because the collective flow re¬ 
sults from AMPT model depend on both the magnitude 
and the distribution of parton scattering cross section. 
Recent analysis m showed that the experiment ellip¬ 
tic flow results can be described by the string melting 
AMPT model with a smaller but more isotropic parton 
scattering cross section. Thus, by carefully tuning the 
parton scattering cross section and its distribution, we 
might achieve a roughly description of all parton species, 
but this is beyond the goal of this paper. 


B. Differential V2{pt)/v3{pt) Ratio from AMPT 
Model 

Figure depicts AMPT model calculations for 
the transverse momentum dependent differen¬ 
tial ratio V 2 {pt)/v 3 {pt) for identified particles 
{Tr^,K^,K'^,(j),p,p,A and A) from 0%-80% central 
Au-|-Au collisions at y'sNN = 200 GeV. The differential 
ratio decreases significantly at low transverse momenta 
{pt < 1.5 GeV/c) and becomes fiat at intermediate 
transverse momenta (1.5 < pr < 3.0 GeV/c). The pt 
dependent differential ratio shows an anti-correlation 
with parton scattering cross section for all particle 
species: the higher the parton scattering cross section, 
the lower the V 2 {pt)/v 3 {pt) ratio. This is in agreement 
with a larger increase of V 3 with increasing parton 
scattering cross section compare to V 2 discussed in Fig. 
and Fig. This effect can also be understood in a 
hydrodynamic frame work. The higher parton scattering 
cross section in the AMPT string melting model is 
equivalent to a lower viscosity in viscous hydrodynam¬ 
ics 001111]. From this study it is known that V 3 is more 
sensitive to viscosity than V 2 m therefore a lower 


viscosity (higher parton scattering cross section) leads 
to a lower V 2 {pt)/V 3 {pt) ratio. 



p.^(GeV/c) 


FIG. 6 . (Color online) AMPT calculation of differential 
V 2 {pt)/v 3 {pt) ratios for identified particles (tt*, A®, (j), 

p, p, A and A) from 0%-80% central Au-|-Au collisions at 
- J snn = 200 GeV. 


C. Integrated V2IV3 Ratio from AMPT Model 

Figure shows the AMPT calculation for the inte¬ 
grated V 2 /V 3 ratio for identified particles from 0%-80% 
central Au-l-Au collisions at ^snn = 39 and 200 GeV. 
The integration region is limited to 1.5 < pt < 2.8 GeV/c 
to avoid the steep decrease in the differential ratio and 
the large error region. Please note, the triangular flow 
used to calculate the V 2 {pt)/v 3 {pt) ratio for 1.5 mb and 
3 mb parton scattering cross section are taken from 0%- 
60% and 0%-70% due to the event plane resolution limits, 
as shown in Fig.j^ The other of ^3 values are taken from 
0%-80% results. 
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FIG. 7. (Color online) Upper panel: AMPT calcula¬ 
tion of the integrated n 2 /v 3 ratio for identified particles 
(tt^, Ag, A and A) from 0%-80% central Au-|-Au 

collisions at ^snn = 39 and 200 GeV. Bottom panel: the 
ratio between integrated W 2 /r >3 of all particles species to that 
of TT^ 


A similar anti-correlation between the integrated f 2 /f 3 
ratio and the parton scattering cross section is observed 
as in the differential V 2 (pt)/v^(pt) ratio: the larger the 
parton scattering cross section, the smaller the integrated 
f 2 /f 3 ratio. In addition, for the same parton scatter¬ 
ing cross section, an energy dependence is observed: the 
lower the energy the larger the integrated V 2 /v^ ratio. 
We further observe that the integrated V 2 /ratio calcu¬ 
lated by string melting AMPT shows no mass (particle 
type) dependence, the ratio between integrated U 2 /U 3 of 
all particle species to that of 7 r+ are consistent with 1 
within 10%. On the other hand the ratios of tt^ from 
default AMPT are lower than the rest of the particle 
species at both ^snn = 39 and 200 GeV. This indicates 
that the quark coalescence mechanism might play an im¬ 
portant role in flattening the U 2 /U 3 ratio. The no mass 
(particle type) dependence of the U 2 /U 3 ratio was pre¬ 
dicted by viscous hydrodynamics, for fast particles (par¬ 
ticles with large px)- The V 2 /v^ ratio is constant and 
does not depend on the mass of the particle [53]. Based 


on the viscous hydrodynamics in Ref. [53|, the collective 
flow depends on dissipative contributions and the Fourier 
coefficient (I^i) of the flow velocity distribution relative 
to the event plane, the dissipative contributions cancelled 
each other in V 2 /vz ratio calculation and only V 2 /V 3 left 
which does not depend on the mass of different particle 
species. Since the parton scattering cross section is re¬ 
lated to the viscosity in hydrodynamics, the magnitude of 
the U 2 /U 3 ratio might be usable to quantify the viscosity 
of the system. 

IV. SUMMARY 

The AMPT calculations of elliptic flow at .^snn = 
11.5, 39 and 200 GeV and triangular flow at y'sNN = 
200 GeV of identified particles (tt^ , , p, p and (p) from 

0%-80% central Au-|-Au collisions were presented. Most 
of results from experiments can be reproduced with a 
parton scattering cross section between 1.5 and 3 mb by 
following the same analysis method as used in the exper¬ 
iments. An anti-correlation between the parton scatter¬ 
ing cross section for both the differential and integrated 
'C 2 /i ’3 ratio is observed. This is due to a different sensi¬ 
tivity of V 2 and U 3 to the parton scattering cross section. 
A mass (particle type) independence of the integrated 
U 2 /W 3 ratio is observed from string melting AMPT, but 
not for the default AMPT. The no mass (particle type) 
dependence of the U 2 /U 3 ratio also can be understood in 
the viscous hydrodynamic framework. Latter one sug¬ 
gests that the magnitude of the ratio might be used to 
quantify the viscosity of the system. 

ACKNOWLEDGMENTS 

We thank Dr. Nu Xu for his great idea and impor¬ 
tant discussions and Dr. Guoliang Ma for his help on 
AMPT model calculations. This work was supported by 
the National Natural Science Foundation of China un¬ 
der grant No.11475070 and U1332125 and the Program 
for Innovation Research of Science in Harbin Institute of 
Technology (PIRS OF HIT B201408). 


[1] S. A. Voloshin, A. M. Poskanzer, and R. Snellings, 
(2008), arXiv:0809.2949 [nucl-ex] 

[2] P. Sorensen, (2009), arXiv:0905.0174 [nud-ex] 

[3] J. Adams et al. (STAR), Nucl.Phys. A757, 102 (2005) 
arXiv:nud-ex/0501009 [nucl-ex] 

[4] K. Adcox et al. (PHENIX), Nucl.Phys. A757, 184 
(2005), arXiv:nucl-ex/0410003 [nucl-ex] 

[5] C. Adler et al. (STAR), Phys.Rev.Lett. 87, 182301 
(2001), arXiv:nucl-ex/0107003 [nucl-ex] 

[6] J. Adams et al. (STAR), Phys.Rev.Lett. 92, 052302 
(2004), arXiv:nucl-ex/0306007 [nucl-ex] 


[7] P. Huovinen, P. Kolb, U. W. Heinz, P. Ruuskanen, 
and S. Voloshin, Phys.Lett. B503, 58 (2001), arXiv:hep- 
ph/0101136 [hep-ph] 

[8] D. Molnar and S. A. Voloshin, Phys.Rev.Lett. 91, 092301 
(2003), arXiv:nucl-th/0302014 [nucl-th] 

[9] J. Adams et al. (STAR), Phys.Rev.Lett. 95, 122301 
(2005), arXiv:nucl-ex/0504022 [nucl-ex] 

[10] B. Abelev et al. (STAR), Phys.Rev.Lett. 99, 112301 
(2007), arXiv:nucl-ex/0703033 [NUGL-EX] 

[11] B. Abelev et al. (STAR), Phys.Rev. C81, 044902 (2010) 
arXiv:1001.5052 [nucl-ex] 










7 


[12] B. Abelev et al. (STAR), Phys.Rev. C75, 054906 (2007) 
arXiv:nud-ex/0701010 [nucl-ex] 

[13] L. Adamczyk et al. (STAR), Phys.Rev.Lett. 110, 142301 
(2013), arXiv:1301.2347 [nucl-ex] 

[14] L. Adamczyk et al. (STAR Collaboration), Phys.Rev. 
C88, 014902 (2013), arXiv:1301.2348 [Unknown] 

[15] B. Alver and G. Roland, Phys.Rev. C81, 054905 (2010) 
arXiv:1003.0194 [nucl-th] 

[16] D. Teaney and L. Yan, Phys.Rev. C83, 064904 (2011) 
arXiv:1010.1876 [nucl-th] 

[17] B. H. Alver, C. Gombeaud, M. Luzum, and J.-Y. Olli- 
trault, Phys.Rev. C82, 034913 (2010), arXiv: 1007.5469 
[nucl-th] 

[18] R. S. Bhalerao, M. Luzum, and J.-Y. Ollitranlt, 

Phys.Rev. C84, 034910 (2011), arXiv:1104.4740 [nucl- 
th] 

[19] L. Adamczyk et al. (STAR), Phys.Rev. C88, 014904 
(2013), arXiv:1301.2187 [nucl-ex] 

[20] X. Sun (STAR Collaboration), Nucl.Phys. A931, 1194 
(2014) 

[21] E. Retinskaya, M. Luzum, and J.-Y. Ollitranlt, 

Phys.Rev. C89, 014902 (2014), arXiv:1311.5339 [nucl- 

th] 

[22] E. Retinskaya, M. Luzum, and J.-Y. Ollitranlt, 

Nucl.Phys. A926, 152 (2014), arXiv:1401.3241 [nucl-th] 

[23] C. Lang and N. Borghini, Eur.Phys. J. C74, 2955 (2014) 
arXiv:1312.7763 [nucl-th] 

[24] A. M. Poskanzer and S. Voloshin, Phys.Rev. C58, 1671 
(1998), arXiv:nucl-ex/9805001 [nucl-ex] 

[25] S. Voloshin and Y. Zhang, Z.Phys. CTO, 665 (1996) 
arXiv:hep-ph/9407282 [hep-ph] 

[26] K. Xiao, N. Li, S. Shi, and F. Liu, J.Phys. G39, 025011 
(2012), arXiv:1111.6213 [nucl-th] 


[27] D. Kharzeev and M. Nardi, Phys.Lett. B507, 121 (2001) 
arXiv:nucl-th/0012025 [nucl-th] 

[28] Z.-W. Lin, C. M. Ko, B.-A. Li, B. Zhang, and S. Pal, 
Phys.Rev. C72, 064901 (2005), arXiv:nucl-th/0411110 
[nucl-th] 

[29] M. Gyulassy and X.-N. Wang, Comput.Phys.Commun. 
83, 307 (1994), arXiv:nucl-th/9502021 [nucl-th] 

[30] X.-N. Wang, Phys.Rev. D43, 104 (1991) 

[31] X.-N. Wang and M. Gyulassy, Phys.Rev. D44, 3501 
(1991) 

[32] B. Zhang, Comput.Phys.Gommun. 109, 193 (1998) 
arXiv:nucl-th/9709009 [nucl-th] 

[33] B. Li and W. Bauer, Phys.Rev. C44, 450 (1991) 

[34] B.-A. Li and G. M. Ko, Phys.Rev. C52, 2037 (1995) 
arXiv:nucl-th/9505016 [nucl-th] 

[35] Z.-w. Lin and G. Ko, Phys.Rev. C65, 034904 (2002) 
arXiv:nucl-th/0108039 [nucl-th] 

[36] L.-W. Chen, C. Ko, and Z.-W. Lin, Phys.Rev. C69, 
031901 (2004), arXiv:nucl-th/0312124 [nucl-th] 

[37] J. Xu and C. M. Ko, Phys.Rev. C84, 014903 (2011) 
arXiv: 1103.5187 [nucl-th] 

[38] S. Adler et al. (PHENIX), Phys.Rev.Lett. 91, 182301 
(2003), arXiv:nucl-ex/0305013 [nucl-ex] 

[39] J.-Y. Ollitranlt, A. M. Poskanzer, and S. A. Voloshin, 
Phys.Rev. C80, 014904 (2009), arXiv:0904.2315 [nucl- 
ex] 

[40] M. Gyulassy, Y. Pang, and B. Zhang, Nucl.Phys. A626, 
999 (1997), arXiv:nucl-th/9709025 [nucl-th] 

[41] B. Zhang, M. Gyulassy, and G. M. Ko, Phys.Lett. B455, 
45 (1999), arXiv:nucl-th/9902016 [nucl-th] 


